Dynamically-Stabilized Pores in Bilayer Membranes 
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Zhelev and Needham have recently created large, quasi-stable pores in artificial lipid bilayer vesicles 
(Biochim. Biophys. Acta 1147 (1993) 89). Initially created by electroporation, the pores remain 
open for up to several seconds before quickly snapping shut. This result is surprising in light of 
the large line tension for holes in bilayer membranes and the rapid time scale for closure of large 
pores. We show how pores can be dynamically stabilized via a new feedback mechanism. We also 
explain quantitatively the observed sudden pore closure as a tangent bifurcation. Finally we show 
how Zhelev and Needham's experiment can be used to measure accurately the pore line tension, 
an important material parameter. For their SOPC/CHOL mixture we obtain a line tension of 
2.6 X 10"''dyn. 

Running title: Dynamically-Stabilized Pores 

Keywords: vesicles, electroporation, lipid, micropipette, line tension 



I. INTRODUCTION 

Lipid bilayer membranes have remarkable physical 
properties. One of the most important among these prop- 
erties is a membrane's resistance to rupture. In the body, 
this resistance is critical to the maintenance of well de- 
fined and properly functioning cells. Indeed, when a cell 
needs to undergo a topological change (as it does during 
cell division, cell fusion, endocytosis, and exocytosis) it 
usually has to make use of specialized machinery which 
carries out the change at the cost of chemical energy. 
This cost is largely determined by the material proper- 
ties of the lipid membranes in question. 

We can quantify a membrane's resistance to rupture 
in terms of a line tension (7), the free energy cost per 
unit length of exposed edge. Edges are disfavored due 
to the high cost of either exposing the hydrophobic lipid 
chains to water, or creating a highly-curved rolled edge 
to hide them. Many authors have devised ingenious in- 
direct measurements of 7 in various lipid systems [Q-D , 
but direct measurement has proven difficult. Among the 
biologically-relevant questions which require such mea- 
surements is the variation of 7 with lipid shape Q . 

Recently Zhelev and Needham have found a new tech- 
nique allowing direct mechanical measurement of the line 
energy |^,^. In this paper we will present a new analy- 
sis of their experimental data. The experiment revealed 
some surprising qualitative phenomena involving pores, 
which we will explain. Briefly (see below), they created 
long-lived quasistable pores about a micron in radius. Af- 
ter persisting for up to several seconds, the pores snapped 
shut in just one video frame. We will quantitatively ex- 
plain the longevity of the pores and their sudden demise, 
fitting several quite different events with a common value 
of 7 and two auxiliary parameters. 

To see why long-lived pores are surprising, consider 
the usual energy of a circular hole in a flat bilayer mem- 



brane Q. This energy can be written as a line tension 
term which is linear in the pore radius minus a surface 
tension term which is quadratic. The energy thus has the 
form 



E{r) — 27r7'7 — Evrr^ 



(1) 



which has only one stable minimum (at r = 0). There 
is a critical radius (r = 7/S) above which the pore 
is unstable to rupture. To cross this critical point, 
the system must surmount a signiflcant energy barrier 
[AE = 7r7^/I]). For typical estimates of the line ten- 
sion (10~^ dyn) thermally-driven rupture thus requires a 
surface tension on the order of one dyn/cm, as observed 
[Q. For lower tensions, any transient pore will reclose 
rapidly, while for larger tensions it will grow rapidly and 
lyse the vesicle; in either case, one does not expect large 
stable pores to exist. 

Nevertheless, Zhelev and Needham found that pores 
with a radius of approximately one micron can remain 
stably open for several seconds. In these experiments, 
giant vesicles are aspirated into the mouth of a mi- 
cropipette where they are held in place by suction. A 
brief electrical field impulse is applied across the vesicle 
by a pair of capacitor plates. As a result of this im- 
pulse, lipid molecules in the membrane rearrange around 
a newly formed pore through a process known as electro- 
poration 

As noted above, the appearance of large stable pores 
is surprising, and yet Zhelev and Needham documented 
over a dozen such events. They proposed that somehow 
these events managed to sit on the unstable equilibrium 
point of Eq. |l| for a long time before suddenly falling 
off. Inferring E and r from the data then let them find 
the line tension from 7 = Er. It seems unlikely, however, 
that the membrane would remain in unstable equilibrium 
for so long. 

In this paper we will find a feedback mechanism that 
dynamically stabilizes pores. The key to the feedback is 
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the relationship between the outflow of solution through 
the pore and the velocity of the vesicle's leading edge as 
it is aspirated into the pipette. The result is a reduced 
surface tension which is a function of both the projected 
length of membrane in the pipette and the pore radius. 
This reduced surface tension yields a new effective energy 
which exhibits a thermodynamically stable pore at finite 
radius. The pore exists for some time before suddenly 
disappearing. From the critical conditions leading to the 
loss of stability, we will be able to produce estimates 
of various parameters in the theory; in particular, we 
will accurately determine the line tension of the bilayer 
membrane. 



II. THE EXPERIMENT 

Figure |l| defines our notation. In the experi- 
ment lipid bilayer vesicles are prepared from stearoy- 
loleoylphosphatidylcohline (SOPC) with 50 mol% choles- 
terol (CHOL) at 24°C. The surrounding solvent is about 
0.5M glucose solution, which effectively prohibits perme- 
ation of water through the membrane by the osmotic 
clamp effect. A micropipette is used to immobilize a 
chosen vesicle using a suction (— p). The suction pres- 
sure is held constant throughout the experiment at a dis- 
tant manometer. Initially a small amount of membrane 
is pulled into the micropipette, leaving a tense spherical 
outer bulb of radius i?init- A square- wave electric field 
pulse is then applied across the vesicle: typical pulses 
produce a field on the order of lOOkV/m for a duration 
of about 150 /zs. The effect of this field is to rearrange 
lipid molecules at one of the vesicle's poles so as to open a 
pore in the membrane. Sometimes no large pore opens; 
in these cases the suction is stepped up and a second 
pulse is apphed. Other times, the electropore is so large 
that the whole vesicle gets sucked rapidly through the 
micropipette and disappears. But occasionally the pore 
stabilizes and the vesicle moves slowly down the pipette 
in a controlled fashion.^] 

What is measured is then the constant applied suction 
(— p) at the manometer, the initial bulb size (i?init), the 
micropipette diameter (dp), and the location {L{t)) of 
the leading edge of the membrane as it advances down 



p=0 



the micropipette. Figure 2 
Other quantities in figure 1 



shows a typical time course, 
such as the bulb radius (i?), 
the pore radius (r), the lubrication layer thickness (/i), 
the surface tensions (S, Eq), and the pressures (— pi, P2) 
are all time dependent and must be inferred from the 
directly measured data. 



'Optical contrast methods make visible the jet emerging 
from the pore, and assure us that the pore does close, rather 
than just getting pulled into the micropipette M. 
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FIG. 1. Geometry of the stabilized pore experiment of 
Needham and Zhelev. 

As the vesicle moves into the pipette the suction im- 
mediately in front of it is reduced due to Poiseuille loss 
along the micropipette. The corrected pressure is given 
by -pi- 



(2) 



Here ry is the viscosity of the sugar solution, L is the 
projected length of membrane in the micropipette (fig- 
ure 1^), and V = L is the velocity of the vesicle's leading 
edge. Eq. || should really be regarded as a definition of 
the effective length (Log) since the micropipette is not re- 
ally a perfect cylinder of constant diameter. Zhelev and 
Needham estimated this parameter as 922 /im by noting 
the velocity at which small beads moved down the tube 
under similar applied pressures. Given the significant dif- 
ference between this experiment and the one in question, 
we will treat L^s as an undetermined experimental pa- 
rameter of fixed value.[] As we shall see later, L^g can be 
determined from the data of figure fusing our theoretical 
approach. 

There is a second velocity-dependent friction term that 
enters into the equation of motion for the vesicle's leading 
edge. This term is due to shear in the lubricating layer 
sandwiched between the membrane and the micropipette 
wall. This frictional force creates a difference between the 
surface tension (Sq) at the leading edge and (S) on the 
exterior bulb: 



(3) 



Here h defines the thickness of the lubrication layer An 
experimental determination of this parameter is some- 
what difficult; Zhelev and Needham estimate it to be on 



^ Left also includes any other constant friction, for example 
membrane drag around the lip of the micropipette. 

■''In reality this parameter will be a function of distance along 
the micropipette and of the projected length itself (R. Bru- 
insma, preprint 1996). It turns out however, that the final 
results are not strongly dependent on h, and so we will treat 
it as a constant. 
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the order of 0.2 — 0.3 fim. As with Leff, we wih determine 
this difficult to measure parameter directly from the data 
in figure ||; our value agrees well with the experimental 
estimate. 
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FIG. 2. Progress of the leading edge of a SOPC/CHOL 
mixed lipid vesicle down the pipette (reproduced from figure 
3of jHj). For this event _Rinit ~ 17.6 /im and p = 353 dyn/cm^. 
The initial velocity (winit) is approximately 35 ^m/s. The solid 
line depicts the theoretical L{t) curve for our stabilized pore 
model. 

These two sources of friction control the speed of the 
vesicle front. As we will show in the following sections, 
they also determine the stable pore size and when and if 
the pore recloses so as to stop the inhalation process. As 
the projected length of membrane inside the micropipette 
grows, so does the amount of friction. We will identify 
a feedback mechanism through which the velocity, and 
hence the surface tension, becomes a function of the pore 
size. This modifies the effective energy function (Eq. |^) 
and in doing so, generates a second stable minimum at a 
finite pore radius (^stable)- At Lcrit the stable minimum 
disappears and the pore snaps shut. 



III. STABILIZATION MECHANISM 

The feedback mechanism mentioned above requires the 
surface tension on the exterior bulb to be a function of 
the projected length of membrane in the pipette and the 
pore size. The first step is to write this tension in terms 
of the pressure inside the bulb (^2) using the Laplace 
formula: 



P2 



2E 



Here R is the radius of the bulb. This formula really 
only applies at equilibrium; fortunately the membrane's 



fast relaxation|j and the fact that R and p2 change rather 
slowly imply that the surface is never far from this ideal. 

The radius of the bulb R is already a simple function of 
L. The necessary relation is obtained from the constraint 
of fixed vesicle areaj^^j We thus have an equation that 
determines the bulb radius as a function of projected 
length (L): 



47ri?2 + nd^L = iirR? 



(5) 



Here i?init is the radius of the bulb at the time the pore 
is created while = dp — h gives the diameter of the 
cylindrical portion of the vesicle inside pipette. Eq. ^ 
assumes the bulb region remains spherical and that the 
area of the leading edge does not change significantly as 
L increases. 

The pressure inside the bulb must still be found in 
terms of L. This is achieved by considering the interface 
at the leading edge of the vesicle. Here we can apply 
Laplace's formula again: 



Pi +P2 



4Eo 



(6) 



Implicit in this equation is the assumption that the pres- 
sure is uniform everywhere inside the vesicle. There will 
in fact be a small pressure gradient near the pipette 
mouth due to convergent flow. This pressure change is on 
the order of rjv/d^ and turns out to be neghgible. Eqs. ^ 
and ^ can now be used to eliminate p2. In addition, we 
can use Eq. |^ to remove So and Eq. ^ to get rid of pi . 
This gives us a new form for the surface tension: 



2R 



2R-d^ 



pdy 
4 



rjv 



8(Leff-i)dv , L 



dl 



(7) 



This equation for the surface tension now has explicit 
dependence on the projected membrane length (both 
through L itself and through the function R{L)). The de- 
pendence on the pore size enters only implicitly through 
the velocity. 

The goal is now to find a suitable expression for the 
velocity to substitute into Eq. |^. Following ||] , this can 
be achieved by considering the flow through the pore. As 
the vesicle moves into the pipette, its volume changes: 



^A rough estimate of the relaxation time is given by 
(4) V « 0.01s; here rj npid ~ lerg s/cm^ is the two- 



dimensional viscosity of the membrane. The pore size (r) 
provides a measure of the disrupted area while the line ten- 
sion (7) gives the restoring force. 

** Strictly speaking, the area is not exactly fixed: thermal 
fiuctuations can always be fiattened out by tension to produce 
a larger projected membrane area. For the tensions appropri- 
ate here, the area change will not be more than a few percent 
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solvent must exit through the open pore. We compute 
the rate of outflow (Q) in two different ways and compare 
to get V. The first way is to write Q as the derivative of 
the total vesicle volume (Figure p: 



d /47r 
"d^ IT^ 



-{2R-d^). (8) 



Note that we have re-expressed all the time derivatives in 
terms of the velocity. The R time derivative was obtained 
by differentiating the area constraint (Eq. ||): 



dR 
'dt 



vdv 
"8R 



(9) 



The other expression for the outflow can be found from 
the pressure forcing the flow. Since the flow is at low 
Reynolds number, Q is just proportional to the pressure 
drop across the pore M : 



3v 



(10) 



Eqs. H and [T^ can now be combined with Eq. |^ to yield 
the desired expression for the velocity of the leading edge 
of the vesicle: 



Sirrjd^ R{2R - d^ 



(11) 



This equation can now be substituted back into Eq. |^ 
to yield the promised form of the surface tension. This 
form has the surface tension as a function of the projected 
length (L) and the pore size (r): 



where 



ai 



a2 



pd^ 2R 



1 + a2r-^ 



and 



4 2i?-dv ' 

16 /8(Leff-iK i 

" h 



(12) 



(13) 



37r(2i? - dv)2dv 



dl 



These equations define the feedback mechanism through 
which the effective surface tension is modified by the fric- 
tion terms outlined in the previous section. The next step 
is to see how this feedback mechanism accounts for the 
experimentally observed behavior — that is, to see how 
it stabilizes pores. 

Our main physical hypothesis is now that at each mo- 
ment the pore adjusts quickly to minimize an effective 
energy similar to Eq. |^, but with the tension replaced by 
the varying quantity just found (Eq. |l^). As long as this 
effective energy has a nontrivial minimum, the pore size 
will track it. This gives the pore size, and hence v = L 
via Eq. |l| in terms of L. We can then solve this ordinary 
differential equation to obtain the time course L(t) and 



compare to the data in Figure This program relies 
on the presence of two different time scales: a slow scale 
for changes of L and R, and a much faster time scale on 
which the pore size r adjusts and the membrane tension 
equilibrates. We are adiabatically eliminating the fast 
variable to obtain a simple dynamics for the slow one. 
Thus our effective pore energy depends on L: 



Ecs{r,L) = 27rr7 



nair 
1 + aiT^ 



(14) 



This new form of the energy can indeed have two minima: 
the trivial one (r = 0) that appeared before and a new 
one at a finite pore size, depending on L: 



dr 



7 







2 (1 + «2rLble)' 



(15) 



As long as L does not exceed some critical value (Lcrit), 
there will be a pore size (rstabio(-^)) that satisfies this 
equation. When the projected length reaches this crit- 
ical value, the second stable minimum disappears and 
the pore collapses (figure ||). At this point the second 
derivative of the energy must also vanish: 



cff 



dr-^ 



= 



(16) 



= 



1 - 7(a2r3j^bic)crit + (a2r3t^bic)crit 



(1 



[air, 



stable 



)crit)3 



The resulting equation is not an equation for the pore 
size, but rather for the product Q;2f stable j ^ pure number. 
The relevant root of equation ( p7| ) yields (a2''stabic)crit ~ 
0.146. This number will be useful in the following section 
where we will use the critical condition to determine the 
three undetermined parameters in the problem: the effec- 
tive pipette length (ioff), the thickness of the lubrication 
layer (ft,), and the line tension (7) of the membrane. 
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FIG. 3. Energy as a function of pore size for supercriti- 
cal, critical, and subcritical lengths of aspirated membrane: 
L = 0,45.8, 90 ^im. (L^wt = 45.8 /xm.) 
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As promised, we have succeeded in reducing the full dy- 
namics of the problem to only one variable, the projected 
length (L). This can be seen by substituting the sur- 
face tension formulae (Eqs. ^ and |l^ along with R{L) 
(Eq. H) and the newly determined pore size {rstahiciL)) 
into the velocity (Eq. |ll|): 



L = v{L) = 



SttT] Vl + Q!2(i)rstable(i)^ 

1 

R{L)d^{2R{L) - d^) ■ 



(17) 



This equation for L as a function of L can now be in- 
tegrated to obtain L{t), then compared directly to the 
experimental results (figure ||). 



IV. EXTRACTION OF PARAMETERS 

In the last section we were able to find a stable pore 
size for each value of the projected length of membrane 
in the pipette. From this pore size we were then able to 
determine the velocity, Eq. This equation is a first 
order nonlinear ordinary differential equation for L. The 
goal of this section is to fit the unknown parameters 7, 
Leff, and h so that the solution L{t) fits the event shown 
in figure |[ We will then use these values to explain the 
other SOPC/CHOL events documented in [||. 

The model that we propose has seven parameters in 
all (table ||). Of these, two are measured directly from 
microscope images, namely the initial bulb radius (i?init) 
and the pipette diameter (dp). The pressure (p) at the 
manometer^ and the viscosity (rj) of the solvent^ are 
also determined experimentally. 

TABLE I. Parameters in the Model 



Name 


Value 


-Rinit, p 


see table ^ 


dp 


7.5 nm 


V 


0.014 erg s/cm^ 


Lc«, h, 7 


fixed by model 



Of the remaining three parameters, Zhelev and Need- 
ham estimated Left and h via auxiliary experiments, then 
deduced 7. As discussed in the second section, the ef- 
fective pipette length was found using a rather different 



^^We determined p from j^: the static tension reported in 
their table I was substituted into their Eq. 4 to yield the 
manometer pressure. 

"'•'''Zhelev and Neeedham used different sugar solutions inside 
and outside the vesicle in order to get visual contrast. For 
simplicity, we just used the average of the two measured vis- 
cosities in our calculations. 



experiment and so in our analysis, it will be determined 
from the data. Our analysis will also yield a value for 
the lubrication thickness {h) that will be in good agree- 
ment with the experimental estimate of Finally we 
will deduce the last remaining parameter in the model: 
the line tension (7). 

Q give the full time course for one event, which 
we will use to find the three undetermined parameters 
listed above. This event is reproduced in Figure ^ and 
appears as the fourth entry in table || below. From 
the experimental L{t) curve we first extract the pro- 
jected length Lcrit = 45.8 /im at pore closure, the ini- 
tial velocity Vmit = 35.5 ^m/s, and the final velocity 
fcrit = 15.7 /zm/s. We will now use these numbers to 
determine the three unknown parameters: the effective 
pipette length (Lcff), the lubrication layer thickness (h), 
and the line tension (7). 

The first step is to require that the pore lose its stabil- 
ity at the observed critical point; the second stable mini- 
mum in the energy must therefore disappear at Lcrit • As 
noted in the last section, the minimum vanishes when 
Eqs. |6| and |l7| are satisfied. To get the pore size just be- 
fore closure (rgtabio.crit) we use Eq. 18 at the critical point 
and require that it reproduce the observed final velocity. 
Recalling Eqs. R4 we then obtain: 



Vcrit 



(«2r,^tablc)crit 



pdl 



where we have defined 



^crit 



) + 



/72 T ■ 

dvh 



(18) 



(19) 



Eq. |l^ can be solved for Ltot- Then from Eq. and the 
value of (a2?'stabie)crit, ouc obtaius the critical pore size: 



^stable, crit 



("2?'stablo)crit 




(20) 



In Eqs. |l| and 1^ the value of R is to be taken from 
Eq. 1^ with L = Lcrit- Using the critical pore size ob- 
tained above, we can now go back to Eq. |l^ to obtain 
a numerical estimate for the product of the line tension 
and a yet to be determined multiplier: 

dp ^ pdpR ( 37rd^(Q2r-3^tabic)crit \ ' 
^dv ~ 4 1^ 16(2i? - dp)Ltot / 

2 - (a2r3^ablo)crit 



(1 



("2r-,3tablc)crit)^' 



' stable icrit 



(21) 



On the right hand side we have made the approximation 
{2R — dv)^/^ « (2i? - dp)^/^ which is accurate to within 
a couple of percent. 

The next step is to determine the parameters Leff and 
h. So far we have found only one combination of these, 
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namely the quantity Ltot introduced in Eq. To sepa- 
rate out the two contributions to Ltot we now require that 
the initial velocity (vinit) come out as observed. From the 
velocity formula (Eq. ^l|) evaluated at L = 0, we can find 
the effective pipette length in terms of the initial pore 
size (rstabio.init) and the initial velocity: 



icff(t'init, ''stable, init) — 



pd. 



37r(2i?i, 



\2j2 



3277i;init 
pdl 3tt{2R 

3277Winit 



1 9Rr^ 

stable, init 



2^2 
P 



128rg^able,init 



(22) 



This form for Leff can then be substituted back into 
Eq. |l^ to obtain a sextic equation for rstabie,init that can 
be solved numerically. The solution so obtained can be 
substituted back into Eq. ^ to yield a numerical esti- 
mate of Leff- Finally the lubrication layer thickness (h) 
can be recovered using Eq. The approximations con- 
sidered in Eqs. |l| and ^ can now be refined using this 
new value of h through a bootstrapping method. 

To finish, we recover the line tension from Eq. ^ by 
multiplying through by the factor (1 — h/dp). Thus all 
three of the required parameters can be determined from 
the single aspiration event depicted in Figure |^: 7 = 
2.6 X 10-*' erg/cm, L^ff = 313 ^lm, and h = 0.48 ^J.m. This 
value of h agrees with Zhelev and Needham's estimate, 
while our Lcs differs considerably from theirs. 

TABLE II. Comparison of experimental 
and theoretical (this paper) closure times for 
SOPC/CHOL vesicles. 



J?init(Mm) 


p (dyn/cm^) 




tth(s) 


Bilayers 


18.3 


399 




2.9 


1 


28.1 


310 


2.6 


3.3 


1 


14.3 


331 






1 


17.6* 


353 


1.7 


1.7 


2 


17.0 


324 


0.33 


0.31 


2 


16.3 


337 


0.33 


0.47 


2 


16.0 


355 


2.0 


1.0 


2 


25.6 


168 


11.0 


11.0 


1 


22.5 


240 


5.5 


5.5 


1 


16.2 


320 


2.5 


4.6 


1 



table compares the experimental time to closure found by 
Zhelev and Needham to the theoretical time determined 
using the model. The same effective pipette length, lu- 
brication thickness and line tension were used for all ten 
events. An effort to optimize the values of the three pa- 
rameters (Leff, 7, and h) based on a least squares fit to 
the critical time data did not produce results significantly 
different from those presented above. 

Although the parameters were fixed by data from a 
single event, they produce reasonable critical times for 
the entire data set. The first event in table Q is clearly 
an exception which we have no explanation for. Perhaps 
there was a large fluctuation in the manometer pressure 
or perhaps the electric field generated a pore so large that 
relaxation to the stable pore size was impossible. Despite 
this anomaly we are confident that we have faithfully de- 
termined the line tension of the SOPC/CHOL membrane 
making up the vesicles in question. 



V. CONCLUSION 

In this paper we have explained the existence of the 
large dynamically-stabilized pores observed by H]. By 
separating the fast timescale on which the membrane re- 
laxation occurs from the slow one associated with the mo- 
tion of the aspirated vesicle down the pipette, we have 
been able to establish a modified pore energy function 
with a stable minimum in the one micron range. In ad- 
dition, we have described a new mechanism by which this 
minimum disappears destabilizing the pore. 

The theory that we have developed permits an accurate 
determination of an important membrane parameter: the 
line tension. From a single event in Zhelev and Need- 
ham's work, we were able to determine this parameter 
and two auxiliary parameters. The values so determined 
were then used to reconstruct all ten of the published 
SOPC/CHOL mixed lipid events: our theoretical post- 
prediction for the critical time at which pore stability is 
lost agreed well with the experimental result for all but 
one event. The agreement supports the value of our the- 
ory as a method for experimentally determining the line 
tension of bilayer membranes. 



* Did not reseal 

^ Event used to determine Lcs, h, and 7 

With the parameters so determined, it is now possi- 
ble to determine the time evolution of the other events 
in 1^ from the initial conditions. For each of the other 
nine SOPC/CHOL events reproduced in table ||, we in- 
tegrated Eq. |l^ for the given initial bulb radius (-Rinit) 
and pressure (p) to obtain the critical time at which the 
pore closed. In some cases a double-bilayer with twice the 
nominal fine tension was needed to fit the data [El. The 
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